Excitations from Filled Landau Levels in Graphene 
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We consider graphene in a strong perpendicular magnetic field at zero temperature with an integral 
number of filled Landau levels and study the dispersion of single particle-hole excitations. We first 
analyze the two-body problem of a single Dirac electron and hole in a magnetic field interacting via 
Coulomb forces. We then turn to the many-body problem, where particle-hole symmetry and the 
existence of two valleys lead to a number of effects peculiar to graphene. We find that the coupling 
together of a large number of low-lying excitations leads to strong many-body corrections, which 
could be observed in inelastic light scattering or optical absorption. We also discuss in detail how 
the appearance of different branches in the exciton dispersion is sensitive to the number of filled 
spin and valley sublevels. 
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I. INTRODUCTION 

Recent experimental progress has allowed the fabrica- 
tion of graphene, a two-dimensional honeycomb lattice 
of carbon atoms that form the basic planar structure in 
graphite Q. Graphene has a host of interesting proper- 
ties descending from its unusual band structure, which 
includes linear spectra around two inequivalent points 
in the Brillouin zone. The eigenfunctions and eigenval- 
ues of the single-particle electron quantum states near 
thesepoints are well-described by the Dirac equation 
|2l IS Li i where the wavefunctions have a spinor structure 
due to the two-point basis of the honeycomb lattice. The 
wavefunctions in the vicinity of each of the Dirac points 
are chiral, leading to a Berry's phase when a particle is 
dragged around either one of them. This has important 
consequences for transport, most notably an absence of 
backscattering that inhibits localization [I| and a finite 
conductivity even when undoped Q . In the presence of a 
magnetic field, the graphene structure has been theoreti- 
cally shown to shift both Shubnikov-de Haas oscillations 
, and the step pattern of the integer quantized Hall ef- 
fect . Both these effects have recently been confirmed 
in experiment [(|. 

All these properties may be understood in terms of 
non-interacting electrons in the graphene bandstructure. 
How do electron-electron interactions impact the prop- 
erties of graphene? To sensibly explore this question re- 
quires an understanding of the low-lying excitations of 
the interacting system. Various parts of these spectra 
may in some cases be measured by inelastic light scat- 
tering [jj or microwave absorption |lCj . In this work, we 
present a study of these excitations in the quantum Hall 
regime. 

For a standard 2DEG, a seminal paper by Kallin and 
Halperin [ll| demonstrated that for densities such that 
the non-interacting spectrum - highly degenerate Landau 
levels with a harmonic oscillator spectrum - has its Fermi 
energy in a gap, the low-lying collective mode spectrum 
may be interpreted in terms of a single particle-hole pair 
excited across the gap. The attractive interaction be- 
tween these particles binds them into an exciton, which 



carries a well-defined momentum P in spite of the mag- 
netic field because the object as a whole is neutral. At 
small P, the exciton dispersion behaves as some form 
of density wave, typically either a magnetoplasmon or a 
spin density wave, depending on the quantum numbers of 
the electron and hole. At larger P, the dispersion crosses 
over smoothly into a form expected for a particle-hole 
pair with separation Pi 2 , with i = yjKc/eB the mag- 
netic length and B the magnetic field. A simple classi- 
cal interpretation of the excitation in this limit is that 
the mutual attraction of the particle and hole at some 
separation precisely cancels the Lorentz force when they 
move together at an appropriate velocity, establishing a 
connection between the momentum of the object and the 
particle-hole separation. 

In graphene this simple picture necessarily breaks 
down, because in the absence of a magnetic field, par- 
ticles and holes move at a constant speed set by the 
linear spectrum around the Dirac points. Nevertheless, 
as we shall see the results for the quantum spectrum of 
these excitations are qualitatively quite similar to those 
of the 2DEG. The main differences arise from the rela- 
tively large number of states within each Landau level - 
there are two spin states and two valley states for each 
Landau state - which allows for sixteen different excita- 
tions when the electron and hole have different Landau 
indices. We shall further see that the large number of lev- 
els leads to stronger many-body corrections to the simple 
two-body model described above than is the case for the 
2DEG. 

In what follows we will examine the exciton spectrum 
of this system in two different ways, first as a two-body 
problem of a positively and negatively charged particle 
each with a graphene spectrum in a magnetic field, and 
then as a many-body problem in which a single particle 
is excited across the Fermi level of a filled Fermi sea. 
At large momentum P we will see these spectra largely 
agree (up to a constant which is undetermined in the first 
approach), while at small P there can be considerable 
difference. 

This article is organized as follows: Section ILTl begins 
with a discussion of non-interacting Dirac particles in a 
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magnetic field. After reviewing the spectrum and wave- 
functions, we turn to the two-body problem, where we 
introduce a canonical transformation to separate the rel- 
ative motion. We solve for the spectrum in the pres- 
ence of the Coulomb interaction, incorporating several 
Landau levels with an eye towards measuring the de- 
gree of Landau-level mixing. In Section IIHI we intro- 
duce the many-body problem with a detailed discussion 
of the low- lying excitations of the non-interacting system, 
for a variety of configurations of the chemical potential. 
The mixing of some of these excitations can be antic- 
ipated from the symmetries of the Hamiltonian. The 
momentum-dependent exciton energies are then calcu- 
lated and catalogued according to spin and pseudospin 
quantum numbers. We conclude by highlighting some 
salient features of the neutral excitations in graphene, 
contrasting them with those of the 2DEG. The calcula- 
tion of matrix elements and dispersions is detailed in the 
Appendix. 



II. TWO-BODY PROBLEM 

The continuum limit of non-interacting graphene is ob- 
tained by adopting a kp approximation |2| about each of 
the two Dirac points, which we label using a pseudospin 
t z with eigenvalues ±1/2, denoted by f|- and -IL For a 
given spin state, the sublattice (A, B) and pseudospin 
(ff, JJ.) degrees of freedom lead to a four-component wave- 
function, which we write in the basis A,^\ B,ij,A,ij,B). 
In addition to the Zeeman splitting, we allow for a much 
smaller pseudospin splitting hw t . The one-body Hamil- 
tonian can be written Hq = H^ n — \gfiBB\s z — hcott z . 
In the absence of a magnetic field, -ffkin 1S given in our 
four-component basis by 



hvp 



( Px +ip y 
p x - ipy 

px - ip y 

V Px + iPy 



(1) 



In a magnetic field, we substitute p — > p — A. In the 
Landau gauge, A = xy/£ 2 , and p y becomes a good quan- 
tum number. For the most part we will use the units 
I = h = Vp = 1, showing these constants explicitly 
when appropriate. We define a dimensionless constant 
(3 = (e 2 / 'el) I '(hv f I '£) to measure the approximate ratio 
of kinetic and potential energies. In these units, the 
Coulomb potential u(r) = e 2 /(er) becomes j3/r. 

For an eigenstate with p y = k, we define the annihila- 
tion operator c = (p x — ix + ik)/y2, and in a magnetic 
field we have 



H kin = V2 
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which describes each valley as a Dirac oscillator. The 
eigenstates of such oscillators are enumerated by an in- 



teger Landau level index which can take any value, pos- 
itive, negative, or zero. We employ the condensed nota- 
tion for single-particle eigenstates 1 = (rti, ki, t\, s\) with 
energy ex, where the respective labels designate the Lan- 
dau level, p v , pseudospin, and spin of the state 1. For a 
given spin state, the eigenfunctions in the Landau level n 
corresponding to t z = +1/2, t z = —1/2 are, respectively, 
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for n = 0, where L y is the length of the system in the 
y direction, and tp n (x) are the usual wavefunctions of a 
simple harmonic oscillator, given in the Appendix. These 
states have energy 

e = V2(hv F /£)sga(n)\n\^ 2 - \gfx B B\s z - hw t t z , (5) 



showing the "v/Ulnl 1 / 2 energy dependence of the Landau 
levels and the particle-hole symmetry characteristic of 
graphene. 

We will frequently use a more compact notation, rep- 
resenting wavefunctions in component form as 



s(n,t,T)ip X (n,t,T)(x - k), 
where r labels the sublattice, and 
sgn(n) 

1 



(6) 



s(n,t,r) 



1 



X(n,t,T) 



; t =ft, t = A 

or t =JJ., r = B 

■ t=t,T = B 

or t =4, t = A 

; t=it,r = A 

or t r = B 

; t=t,r = B 

or t=ij.,r = A 



(7) 



(8) 



For the two-body problem, we neglect Zeeman and 
pseudospin splitting and consider only the kinetic en- 
ergy of the electron (1) and hole (2), described by the 
Hamiltonian 



H = n£l(A)+Hgl(-A)-u(r 



(2), 



.(1) 



,(2) 



(9) 
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where the sign of the vector potential is reversed in the 
hole term to account for its opposite charge. Without loss 
of generality, we may assume that the electron and hole 
are in the § valley. One may obtain other cases by inter- 
changing the sublattices of one or both particles, which 
does not change the spectrum since the kinetic energy is 
degenerate in the two valleys and the Coulomb interac- 
tion does not scatter between valley states in our con- 
tinuum approximation |12| . With the pseudospin fixed, 
we need only consider the upper left block of H^ in in Eq. 
^ which can be conveniently expressed using the Pauli 
matrices o~i and the 2x2 identity matrix /. On the other 
hand, 4x4 matrices are required to describe the action 
of the Hamiltonian on the sublattice degrees of freedom 
of the two-body system, and we express these matrices 
using a tensor product of two 2x2 matrices in which 
the left and right slots correspond to the action on the 
electron and hole spinors, respectively. We define ccntcr- 
of-mass and relative coordinates as = P/2 + p, 

p (2) = p/2 _ P: r (i) = r + r /2, r (2) = R r /2. We 

now have 

H = P x m+ + (P y - x)m+ (10) 
+2p x m x ~ + 2(p y - X)m~y - u(r) I 

where m x = \{a x ®I±I®a x ), = —~{a y ®I±I®Gy), 
and I = I <g> I . To decouple the center-of-mass motion, 
we introduce the transformation U = e zXy , for which 
U f P x U = P x +y and U*p y U = p y + X. The result is 

U ] HU = 2p ■ m + (P - z x r) • m+ - u(r) I (11) 

It is now apparent that P commutes with the trans- 
formed Hamiltonian and can be replaced by a center- 
of-mass momentum quantum number. For a given P, 
we may study the interacting problem in the relative co- 
ordinate since Wu{r)U = u{r). The operator UPW 
commuting with the original Hamiltonian is identical to 
the two-body momentum operator Q found in Ref. ITU 
When P — 0, the system possesses an additional sym- 
metry consisting of a simultaneous relative spatial coor- 
dinate and SU(2) sublattice rotation. We will return to 
this later. 

We may now shift r — > r — z x _P, moving the in- 
dependence to the potential energy. The kinetic en- 
ergy is now 2p ■ m — (z x r) • m + , which describes 
a two-dimensional Dirac oscillator. To see this, de- 
fine the harmonic lowering operators a — p x — ix/2, 
b = p y — iy/2. We can also take the combination 
c± = (a ± ib)/y/2, which describes independent "+" 
and "-" oscillators. Using the four-component ba- 
sis {a<UaV\AWbP\bWaP\bMbW) for the sublat- 
tices of the electron and hole, 



or alternatively, in our tensor product notation, 



2p m — (z x r) 
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,(12) 
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(13) 



It is now clear that the two-dimensional Dirac oscillator 
is characterized by two quantum numbers n + and n_, 
which are the Landau level index of the electron and 
hole, respectively. 

The eigenstates of the kinetic energy can be written in 
terms of the two-dimensional harmonic oscillator wave- 
functions $ nii „ 2 (r), defined in the Appendix. For an 
electron in Landau level n + and hole in level n_, the 
four-component wavefunctions for the relative coordinate 



( s+s_*[ n+ |_ li[n _|_i(r) \ 
s+$| n+ l_i,l„_l(r) 



V 



ln+Wn-l-1 



(r) 



(14) 



where s± = sgn(n±). The corresponding kinetic energy 
is V2 [sgn(ri + )|n + | 1 / 2 - sg^n^n^ 1 / 2 ]. 

We pause to consider the virtues of the canonical trans- 
formation. It maps the space of momentum states of an 
electron and hole, each of which carries a one-dimensional 
momentum in our gauge, to a two-dimensional space of 
momentum states labelled by P. The relative motion is 
described by a two-dimensional Dirac oscillator moving 
in a _P-dependent potential, whose two non-interacting 
quantum numbers are the Landau level indices of the 
electron and hole. It is worth noting that the separation 
of the center-of-mass coordinates relies on the uniformity 
of the magnetic field and on the overall neutrality of the 
system. When these conditions are met, such a transfor- 
mation is always possible, although it will vary depending 
on the choice of gauge. It has an essential two-body na- 
ture and cannot be reduced to a choice of single-particle 
basis. 

Furthermore, the transformation is not peculiar to 
graphene and is equally applicable to the 2DEG, where 
one obtains a two-dimensional harmonic oscillator. In 
this case, the operators c+ and c_ operate on indepen- 
dent clockwise and counter-clockwise modes. As a result, 
the angular momentum of the oscillator corresponds to 
the difference between the Landau level indices of the 
electron and hole; indeed, ( E > n+ .n_('") has azimuthal de- 
pendence e - l ^(™+-™-) j where cf> is the azimuthal angle of 
r. This motivates the nomenclature l z = n + — n_ for the 
difference between the Landau level indices of the elec- 
tron and hole. When P = 0, the effective potential for 
the oscillator u(r — z x P) becomes rotationally symmet- 
ric, and matrix elements between two-body states with 
different values of l z will vanish. The conservation of l z 
at P = can be viewed as a conservation law for the 
internal angular momentum of the exciton. 

An analogous conservation law arises in graphene. 
However, l z is not conserved, since in an eigenstate of 
the kinetic energy, the different spinor components can 



4 



have different azimuthal dependences. We find instead 
that when P = 0, the quantity |n+| — |n_| is conserved 
by interactions. Physically, this is due to the aforemen- 
tioned P = symmetry involving combined coordinate 
and sublattice rotations. This is the grapheme analog of 
l z conservation in the 2DEG. 

To summarize, the dispersion of the two-body prob- 
lem is obtained by diagonalizing the relative coordinate 
Hamiltonian, whose matrix elements h are given by 

f l n + ,n_;n' + ,n'_(P) = (15) 
* n+ ,„;*n_,n'_V^(« + |n+| 1 / a - S >-I 1/2 ) 

— [ A /2]' 5 "+ 0+ ' 5 ' l -' 0+<5 "+- 0+ ' 5 "'-. ~ 4 . f 

S + S-S' + s'_U\ n+ \-i 1 \ n _\-i ; \ n ' + \-l ! \ n '_\-i 
+ S+S+W| n+ |-l,|n_|;|ri' + |-l,|ri'_| 
+ S-S_1t|n + |,|n_|-l;|n^_|,|n'_|-l 

+ W|„ + |,| n _|;| n ^|,|„'_|] 

where 

u n 1 ,n 2 ;n 3 ,n4 {P ) = (16) 

/ dru(r- zx P)$* 3 ,„ 4 (r)$ 

The location of the chemical potential will determine 
the possible Landau level indices for the electron and 
hole. The splitting of each level into 4 spin and pseu- 
dospin sublevels enlarges the realm of possibilities. We 
use the notation 7n for the number of filled sublevels in 
Landau level n, where n refers to the highest Landau 
level with filled sublevels. Thus, the undoped system 
with chemical potential \x = corresponds to 70 = 2, 
whereas placing the chemical potential between Landau 
levels n and n + 1 is denoted as 7 „ — 4. We consider the 
two-body problem both for j n = 4 and 7 „ < 4, neglect- 
ing the contribution of Zeeman and valley splitting to the 
non-interacting energy. The cases 7„ = 1, 2, 3, which are 
essentially identical in this two-body analysis, raise the 
possibility that both the particle and hole are in level n, 
giving rise to intra- Landau level excitons carrying spin 
or pseudospin. 

Notably, the typical ratio of kinetic and potential en- 
ergies (3 is of order unity and independent of the strength 
of the magnetic field. We would like to establish whether 
Landau level mixing is weak in this situation, partic- 
ularly in order to justify approximations used in the 
many-body treatment of the next section. We have 
(3 = (e 2 /e£)/(hv F /£) = a(c/v F e), where a w 1/137 is the 
fine-structure constant. Experiments both in transport 
PI 0| and far-infrared absorption 0] find c/vf ~ 300. 
Alicea and Fisher ^| estimate e = 5 for graphene grown 
on a SiC>2 substrate, which accounts both for screen- 
ing currents in the substrate as well as weak intrinsic 
screening in graphene computed in the RPA 14] . The 
screening should be somewhat better on a SiC substrate. 



We assume the more pessimistic value e = 3, leading to 
(3 = 1/1.37 w 0.73. 

We study the cases j n = 4 and 7„ < 4 for n — 0, 1, 2, 3, 
computing the dispersion of the first five two-body states 
up to P£ —12 using 4 electron and 4 hole levels. The 
ket I n + ,n_ ) denotes a non- interacting state with the 
electron in level n+ and the hole in level n_. We will 
see that, even though (3 is of order unity, Landau level 
mixing plays a small role for these excitations. Figure 
H shows the case 7n = 4. The dispersions approach the 
non-interacting energies of the states | n + , n_ ) , shown as 
dotted lines, as P — > 00. An analysis of Eq. 1161 shows 
that the limit is approached as 1/P. For n — 0, the non- 
interacting degeneracy in the third excited state (between 
I 4,0) and | 1,-1)) is lifted at low P by interactions. 
The dispersions computed from a single state | n+, n_ ), 
without Landau level mixing, are shown as dashed lines. 
Although (3^1, the effect of mixing is relatively weak for 
many of the exciton states. The first and second excited 
states for n = 2 and n = 3 are notable exceptions. This 
can be understood by noting that the non-interacting 
states I n + 1, n — 1 ) and | n + 2, n ) can admix strongly 
at small P since they have the same nearly-conserved 

quantum number |n+| — |ra ] . For the ground state, the 

main qualitative effect of Landau level mixing in all four 
cases seems to be the deepening of peaks and dips in the 
dispersion. 

The same observations apply to the case 7„ < 4, shown 
in Figure Here, there are non- interacting degeneracies 
in the excited states for n = 0, and in the third excited 
state for n = 1. Note that in the former case, this de- 
generacy is lifted at low P only after the inclusion of 
Landau level mixing. The degeneracy is not broken for 
P = 0, at which \n + \ — |n_| is exactly conserved. Due 
to the small (possibly zero) kinetic energies, some exci- 
ton energies become negative. This situation is rectified 
in the many-body calculation of the next section when 
self-energy effects are included. 

We pursue Landau level mixing further in Figures |3| 
and 21 which show the square-modulus of the overlap 
between the ground state (denoted | )) and higher non- 
interacting states, given by 

1-|(0 I n + l,n)| 2 ; 7n = 4 
1-|(0 I n,n)\ 2 ; ln < 4. [U > 

The level of mixing is 15% or less for all cases, vanishing 
rapidly in the limit P — > 00. Maximum mixing occurs 
at higher momenta as more Landau levels are filled. We 
caution, however, that mixing can be much larger for 
even modest increases in (3. For example, we find mixing 
of up to 70% at intermediate momenta for (3 = 2.2. 

The peaks and valleys apparent in the exciton disper- 
sions of Figuresnand[21are of particular interest, because 
they lead to van Hove singularities in the collective mode 
density of states. Such singularities are in principle de- 
tectable in inelastic light scattering and, via coupling 
through disorder, microwave absorption |l5j |. Thus it is 
useful to know when and where such singularites might 




(c) 72 = 4 

FIG. 1: Dispersion of the first few levels for (3 — 0.73 when 7, 



(d) 73 = 4 

= 4. Dashed curves neglect Landau level mixing. 




(c) 72 < 4 

FIG. 2: Dispersion of the first few levels for (3 



(d) 73 < 4 

= 0.73 when 7„ < 4. Dashed curves neglect Landau level mixing. 
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FIG. 3: Amount of Landau level mixing at momentum P, 
as denned by Eq. 1171 for — 0.73. Results shown for 70 = 
4 (solid), 71 = 4 (dashed), 72 = 4 (dot-dashed), 73 = 4 
(dotted). 
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FIG. 4: Amount of Landau level mixing at momentum P, 
as defined by Eq. 1171 for — 0.73. Results shown for 70 < 
4 (solid), 71 < 4 (dashed), 72 < 4 (dot-dashed), 73 < 4 
(dotted). 

be found. To compute them more accurately, we now 
turn to the many-body formulation of this problem. 



III. MANY-BODY PROBLEM 

We study the many-body exciton dispersion by diag- 
onalizing the many-body interacting Hamiltonian on a 
subspace of trial wavefunctions, which consist of a sin- 
gle electron excited from the Fermi sea. In the 2DEG, 
this procedure is known to give results for the exciton 
dispersion that coincide with those obtained from the 
poles of response functions computed in a generalized 
random phase approximation (RPA) 0| . Indeed, these 
wavefunctions are exact in the non-interacting case, and 
should thus provide a good basis for the weakly inter- 
acting regime. On the other hand, the validity of these 
wavefunctions may persist for stronger interactions if the 
Landau level mixing is weak. Experience has shown that 



in the 2DEG, the degree of Landau level mixing is con- 
siderably less than one would naively anticipate based on 
the ratio of potential and kinetic energies. This mixing 
is known to be small 01 even when this ratio is of order 
unity. 

In the case of graphene, this notion is supported by 
the results of the previous section. In the two-body cal- 
culation, where we take (3 = 0.73, we find that in the 
lowest exciton state, Landau level mixing occurs at a low 
level, causing the inflections in the dispersion to become 
somewhat more pronounced, but otherwise having little 
qualitative effect. From this we conclude that reasonably 
accurate results for the lowest exciton states may be ob- 
tained using only those electron and hole Landau levels 
nearest the Fermi energy. Thus, we focus our attention 
on exciton states in which the Landau level indices of the 
electron and hole differ by at most l z = 1. 

As in the previous section, we make the approximation 
that the Coulomb interaction has SU(2) symmetry with 
respect to the two sublattices. In this approximation, 
the Coulomb interaction conserves pseudospin in addi- 
tion to the usual spin conservation. As a result, exciton 
states mixed by the Coulomb interaction have the same 
net s z and t z and consequently have the same net contri- 
bution to the non-interacting energy from Zeeman and 
pseudospin splitting. 

We restrict our consideration of mixing to those ex- 
citon states with the same non-interacting energy. In 
fact, such degeneracy is almost always due to the sev- 
eral possible spin and pseudospin configurations of the 
individual electron and hole when the net s z or t z is 
zero. Due to the nonuniform spacing of Landau levels 
in graphene, two electron-hole pairs with different Lan- 
dau level indices will not have the same kinetic energy 
unless one of the pairs has l z > 1. The only exception 
occurs when 70 < 4, in which case the non-interacting 
states I 1, ) and | 0, — 1 ) have the same kinetic energy. 
Under these considerations, we concentrate on the cases 
70 < 4, 70 = 4, 71 = 4, and 72 = 4. In particular, 
the first case offers the possibility of Landau level mixing 
for l z — 1 and intra-level spin and pseudospin waves for 

h = 0. 

At this level of approximation, the dispersion relation 
E(P) for an exciton can be written 

E(P) = E kin - \gfi B B\s z ~hu t t z + AE U {P), (18) 

where Skin, s z , and t z are, respectively, the combined 
kinetic energy, spin, and pseudospin of the electron and 
hole. The superscript v provides an arbitrarily assigned 
index for the several branches of the additional contribu- 
tion AE"(P) due to interactions. 

In our many-body approach, AE has three contribu- 
tions. The first is the direct Coulomb interaction between 
the excited electron and hole, whose matrix elements are 
identical to those found in the two-body approach of the 
previous section. This "excitonic effect" can be obtained 
in a diagrammatic expansion by including vertex correc- 
tions |ll|. In our continuum approximation, this contri- 
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FIG. 5: Exciton states for the case 7 = 4, with definite spin 
and pseudospin denoted (s z ,t z ). The dashed box indicates 
states admixed by exchange interactions. Direct Coulomb 
interactions do not admix these states. 



bution only scatters between exciton states if the elec- 
tron and hole individually have the same s z and t z in 
the initial and final state. The second is the exchange 
term or "depolarization effect," which is the only contri- 
bution present in RPA. Physically, it represents the an- 
nihilation of the electron-hole pair at one point in space 
and its recreation at another point. This process only 
affects excitons with vanishing net s z and t z . Finally, 
each single-particle sublevel has an associated self-energy 
which arises from exchange interactions with those occu- 
pied sublevels in the Fermi sea having the same spin and 
pseudospin. This effect provides a P-independent shift 
to the exciton dispersion, guaranteeing that its energy is 
always non-negative. It also pushes the P — > oo limit 
above the non-interacting energy. 

The situation is illustrated in Figure [5] for the simplest 
case 7„ = 4. There are 16 exciton states with minimal 
kinetic energy which are differentiated by the spin and 
pseudospin sublevels of the electron and hole. Direct 
interactions only scatter electrons (holes) into electron 
(hole) sublevels with the same s z and t z , and thus do not 
cause mixing. Exchange interactions scatter among the 
4 states with s z — t z = 0. 

The case 7 < 4, l z = 0, is shown in Figure El The 
lowest exciton states actually have zero kinetic energy; 
the excitation is a pure spin or pseudospin wave within 
the Landau level. Neither the direct nor exchange inter- 
actions mix the states. The dispersions are identical to 
that of the two-body case up to a constant due to the 
exchange self-energy correction as well as Zeeman and 
pseudospin splitting. 

The most interesting scenario is illustrated in Figures 
I7l9l for the cases 70 = 1,2, 3. Since the levels n = — 1, 0, 1 
are equally spaced, there are excitons with l z — 1 whose 
kinetic energies are degenerate and may be split by Lan- 
dau level mixing. This is caused by the direct Coulomb 
interaction, which admixes states connected by a wavy 
line. As in the 7« = 4 case, the s z = t z = excitons 
are mixed by the exchange interaction. Note that the 
number of exciton states with given values of s z and t z 
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FIG. 6: Exciton states for the case l z — 0, 7 < 4, labeled by 
spin and pseudospin as (s z ,t z ). 
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FIG. 7: Exciton states for the case 70 = 1, l g = 1, with 
definite spin and pseudospin denoted (s z ,t z ). The dashed 
box indicates states admixed by exchange interactions, wavy 
lines show states coupled by the direct Coulomb interaction. 



is closely tied to the number of filled sublevels. 

We now derive the three many-body contributions by 
diagonalizing the Hamilonian on the subspace of trial 
wavefunctions. The matrix elements between the non- 
interacting exciton states are 

M 1A1 , 2 - = ( O I (a\,a 2i y Ha\a 2 \ Q ) (19) 
-(0|iT|n)<W 2 o,, 



where | Q ) is the filled Fermi sea. Diagonalizing H on 
the subspace of trial wavefunctions may be viewed as a 
two-body problem with Hamiltonian M. The many-body 
Hamiltonian has the form 



H 



E 



- E 

1,2,3,4 



Vi, 2,3, 4a\ala 3 a 4 , (20) 



where Vi ; 2,3,4 = ^2,1,4,3 due to the inversion symmetry of 
the Coulomb interaction. This leads to M = M° + M d + 
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FIG. 8: Exciton states for the case 70 = 2, l z — 1, with 
definite spin and pseudo-spin denoted (s z ,t z ). The dashed 
box indicates states admixed by exchange interactions, wavy 
lines show states coupled by the direct Coulomb interaction. 



FIG. 9: Exciton states for the case 70 = 3, l z = 1, with 
definite spin and pseudospin denoted (s z ,t z ). The dashed 
box indicates states admixed by exchange interactions, wavy 
lines show states coupled by the direct Coulomb interaction. 



M 2 



M n where 



2;1',2' 



Ml 



M d . , 

iw l,2;l',2' 



2;V, 2' 



iW l,2;l',2' 



Si,y5 2 ,2>(ei-£2) (21) 

-Vl',2,2',1 (22) 

V V ,2,1,2> (23) 

^ 6(M- £3)^1,1^3,2,3,2' (24) 



— 62,2' Vi',3,1,3]. 

These terms represent the noninteracting, direct, ex- 



change, and exchange self-energy contributions, respec- 
tively. We assume that the direct interaction with the 
Fermi sea is cancelled by a neutralizing background 
charge. The calculation of these matrix elements is 
greatly simplified by separating the relative and center- 
of-mass motion of the electron and hole by means of 
the canonical transformation developed for the two-body 
problem. The computation of matrix elements and 
AE(P) is discussed extensively in the Appendix. 

For a given chemical potential, interactions give rise 
to several branches AE U {P). Since there are fewer dis- 
tinct branches than there are states (16 states in the case 
l z = 1, 4 for l z = 0,) some of the exciton states remain 
degenerate. Some of this degeneracy may remain even 
after pseudospin and Zeeman splitting are taken into ac- 
count, which is readily done by inspection of the s z and 
t z quantum numbers in Tabled Based on an analysis of 
Figures 15191 and of the matrix elements, we can deduce 
the different spin and pseudospin quantum numbers car- 
ried by excitations with the same AE(P), as is shown 
in Table [I] Analytic expressions for these branches are 
given in the Appendix in Table ITT1 

Along with s z and t z , we may also consider total spin 
s 2 = S(S + 1) and total pseudospin i 2 = T{T + 1). The 
Fermi sea must be in a spin (pseudospin) singlet state in 
order for S (T) to be a meaningful quantum number for 
the excitons, as has been pointed out previously for the 
2DEG 11]. Consequently, excitons have well-defined S 
only when j n = 4 and well-defined T only when j n = 2, 4. 

From symmetry considerations, we can predict an ef- 
fect of intervalley scattering. Such scattering will not 
pick out a preferred direction in pseudospin space, but 
may lift any degeneracy between pseudospin singlet and 
triplet states. Based on Table U this effect should be 
easiest to discern in the case 70 = 2. In this case, there 
are spin-carrying excitations (both for the intra-level ex- 
citon as well as branches AE 3 and AE 4 of the l z = 1 
exciton) having four-fold degeneracy purely due to pseu- 
dospin (provided fno t = 0.) An estimate on the strength 
of intervalley scattering could be gleaned from any split- 
ting observed between the T = and T = 1 states. 

Figures HUI- 1131 plot the various branches AE V {P) for 
the different cases. The l z — case (intra-level exciton, 
Figure HOfl is the simplest; it has a single branch AE° 
given by the two-body result plus a constant exchange 
self-energy correction. This correction renders the mode 
gapless (if we ignore Zeeman and pseudospin contribu- 
tions to the non-interacting energy,) and we recognize it 
as the quadratic Goldstone mode corresponding to bro- 
ken spin and pseudospin symmetries. Next in complexity 
are the cases 71 — 4 and 72 — 4, which include an ex- 
change interaction among the s z = t z = states. This 
interaction splits off the totally symmetric combination 
with s z = t z — and S = T = 1 in branch AE 1 , and 
the other 15 states are in branch AE 2 . The former has a 
nonmonotonic character (particularly for larger n,) while 
the latter is essentially monotonia 

Figure IT31 shows the most interesting case n = 0. The 
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TABLE I: Quantum numbers of the various exciton branches 
for different numbers of filled sublevels. Dash indicates the 
quantum number is not well-defined. For j n = 4, branch AE 2 
contains all the spin and pseudospin states not contained in 
branch AE 1 . In the 70 = 2, l z — case, all 4 pseudospin 
states are present. 

case 70 = 4 is similar in all respects to the other cases 
7„ = 4, having 15 states in branch AE 2 and a single state 
pushed up by exchange interactions into branch AE 1 . 
When less than four sublevels are filled, two additional 
branches appear due to the possibility of Landau level 
mixing. In general, we expect the cases 7„ = 1,2,3 to 
have different P dependences for AE V (P) since the ex- 
change self-energy parameters (which generally do not 
simply add a constant shift to the energy) depend on the 
number of filled sublevels. However, for n = 0, there is 
a degeneracy between exchange constants which is pro- 
tected by particle-hole symmetry (this point is elaborated 
in the Appendix,) and the branches AE"(P) are identical 
for these three cases. This feature is peculiar to graphene 
has no analog in the 2DEG. 

On the other hand, the number and type of excita- 
tions in each branch does depend on the particular value 
of 70, as can be seen in Table |T| Note that if we were 
to add the constant Zeeman contribution to Figures ITUl 
1131 we would see branch AE 2 separate into three evenly 
spaced levels for the cases 70 = 1,3 and 7„ = 4, but not 
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FIG. 10: Energy shift AE for the intra- Landau level exci- 
tons in the case 70 < 4. 
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FIG. 11: Energy shifts AE 1 (solid) and AE 2 (dashed) for 
71 = 4. 



70 = 2. For 70 = 1,3 we would also see branches AE°, 
AE 3 , and AE 4 split into two levels. These effects can 
be deduced from the s z values listed in Table [I] Simi- 
lar considerations apply to the much smaller pseudospin 
splitting. Note that even after the inclusion of these ef- 
fects, the dispersions of the 70 = 1, 3 cases are identical, 
due physically to particle-hole symmetry. 

We remark that the only P-dependent contribution to 
the dispersion not captured in the two-body analysis ap- 
pears in branch AE , present in the l z = 1 excitons. 
Here, exchange interactions lead to dramatic differences, 
for example as is evident in Figure 1121 The very high 
peaks apparent there not only do not appear in the two- 
body formulation, they are quantitatively much larger 
than analogous many-body corrections for the 2DEG 
[Tl| . This relatively large effect may be traced to the un- 
usually large number of pair excitations coupled together 
in the graphene system. It would be most interesting 
if the peak/dip structure evident in Figures ITT1 and fP2l 
could be observed in inelastic light scattering or optical 
absorption. [Toj 
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FIG. 12: Energy shifts AE 1 (solid) and AE 2 (dashed) for 
72 = 4. 
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FIG. 13: Energy shifts for inter- Landau level excitons in the 
case n = 0. When 70 = 4, only AE 1 and AE 2 (solid) are 
present. For 70 < 4, the additional branches AE 3 and AE 4 
(dashed) appear. 



IV. CONCLUSIONS 

In this article, we have analyzed the excitations from 
filled Landau levels in graphene using two approaches. 
In the first, we tackle the two-body problem for the ex- 
cited particle and hole in a magnetic field. We introduce 
a canonical transformation which separates the relative 
motion and allows us to identify the symmetries and con- 
servation laws for the system. We then compute the dis- 
persion of excitons as a function of P, restricting the 
available electron and hole levels based on the value of the 
chemical potential. By including higher Landau levels in 
the calculation, we establish that Landau level mixing is 
weak for the lowest exciton states even though the esti- 
mated ratio P of potential and kinetic energies is of order 
1. We caution, however, that mixing can become impor- 
tant for even modest increases above the value (3 = 0.73 
used here. 

Armed with this insight, we turn to the many-body 
problem, where we consider the implications of the spin 



and pseudospin sublevels. We find that the large num- 
ber of degenerate states coupled together by exchange 
interactions gives rise to an excitation with enhanced 
many-body corrections and strong peak/dip features in 
its dispersion, which could be observed experimentally. 
A careful analysis of the spin and valley quantum num- 
bers shows the sensitivity of the spectrum to the number 
of filled sublevels 70 in Landau level n — 0. In particular, 
when 70 < 4, there are new branches in the dispersion 
relation (due both to Landau level mixing and to intra- 
level excitations) which do not appear when 70 = 4. We 
also find that the case 70 = 2 is distinguished from the 
other cases by the absence of any Zeeman splitting of the 
low-lying exciton states. The cases 70 = 1,3 arc identi- 
cal due to particle-hole symmetry effects which have no 
analog in the 2DEG. 
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APPENDIX A 

The oscillator wavefunctions for the one-particle basis 
defined in the text are 

M*) = (V^2 n n\y 1/2 H n {x) e"* 2 / 2 , (Al) 

where H n are Hcrmitc polynomials. The two- 
dimensional oscillator wavefunctions are obtained by the 
repeated application of the raising operators c± to the 
ground state. Writing $ 0tQ (z,z) — [2Tz)~ l / 2 e~ zz ^ i with 
z = x + iy, we have $ ni .n 2 °c [c+] ni [cl]" 2 $o,o ; or explic- 
itly 

* ni ,„ 2 (r) = (2 7 r)- 1 / 2 (2 ni+ " 2 n 1 ! n 2 \ )" 1/2 (A2) 
{-2d z + z/2) ni {2d- z - z/2)™ 2 e- zS ^ 

= (2 7 r)- 1 / 2 2-l'-l/ 2 -^— e-^sgn^) 1 ' 

r \i*\Ll°}{r 2 /2)e- r2 l\ 

where L denotes Laguerre polynomials, z = x + iy, l z = 
m — n,2, n_ = min(ni,7i2), z/\z\ = e 1 ^ and sgn(Z z )'* — > 
1 for l z = 0. The one- and two-dimensional oscillator 
wavefunctions are connected by the relation 

*n 1 ,n 2 (r) = (-ir+"» / ^=e iX y (A3) 

Vni {X-x/2)<p n2 (X + x/2). 

The separation of the relative coordinate of the two- 
body state can be achieved with an appropriate choice 
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of the two-body basis. We choose the basis in which elements are 
the electron and hole lie in Landau levels n\ and n 2 , 
respectively, and the pair has definite momentum P. For 
the 2DEG, the one-body eigenstates are defined 



^ D fc EG (r) = L-V*e ik y<p n (x - k). (A4) 



The two-body wavcfunctions are formed using a special 
linear combination of one-body states for the electron 
and hole: 



«„ 2 (n,r 2 )^ 



dk 

* 2DEG 



(A5) 



? iXy e ip - R $ nun2 (r + zxP), 



where r and R are the relative and center-of-mass coordi- 
nates corresponding to r% and r%. This basis implements 
the canonical transformation developed for the two-body 
problem. The second equality can be established using 
Eq. HI 

The same transformation connects the one- and two- 
body bases for graphene, which we review here in com- 
ponent form. For the one-body states, 



M 



1,2:1'. 2' 



Mf, 



m l,2;l',2' 
O 



Tl,T 2 ^ 

u(n - r 2 )[^P'.y 1 y 2 ,t[t' 2 (ri,r 2 )}* Ti T2 

[*P:n 1 ,n 2 :t 1 t 2 (ri,r 2 )] ri , T2 

= S SuS2 5 s ^ s ' 2 5t u t 2 St' 1 ,t! 2 '^2 / dridr 2 
u{r x -r 2 )[*P';„' l ,„' 2 ; t ' 1 t' 2 (ri,ri)]; i!ri 

r 2 ,T 2 

X SI c n 

= 0i,i/m 2) 2' - <J 2 ,2''m ljl ,, 

= /] ®L» - £ (^) s l> h)] SaxJ, &t u f, 
n 

ri,r 2 



W n u{x,y)] T = L-^^W~2\ 



(A6) 



s(n,t,r)¥> A ( n) t )T )(a; - k), 



and for the two- body states, 
[*f ;ni ,„ 2 ; tl ,t 2 (ri,r2)] Tl , T2 

/,*2 



f/A- 



iP^k 



w 2 



(r 2 )]; 2 [v 



Si 



/2tt 



(A7) 



e^e lP - fl [v / 2] 5 "i' 0+ ^- - 2 s(ni,ti,T 1 ) S (n 2 ,i 2 ,r 2 ) 
• *A(T. l ,t 1 ,n),A(» J ,( 2 ,7i 1 )('' + z x P), (A8) 



where t\ and i 2 denote the pseudospin index of the elec- 
tron and hole sublevels. 

The separation of the center-of-mass and relative co- 
ordinates in the latter expression is extremely useful in 
the calculation of matrix elements of M. The matrix 
elements Mi )2; i' j2 < represent the scattering between ini- 
tial and final states for the electron- hole pair. Let n\, 
n 2 , n'i, n' 2 denote the Landau level indices of the initial 
electron, initial hole, final electron, and final hole states, 
respectively. Similar conventions denote the initial and 
final momentum P of pair as well as the spin s and pseu- 
dospin t of the individual electron and hole. The matrix 



where u(r) = j3Jr is the Coulomb potential. The inte- 
grals are straightforward after changing to relative and 
center-of-mass coordinates. For the direct (M d ) and ex- 
change (M x ) contributions, the use of Eq. IA8I leads 
immediately to expressions involving $, whereas for the 
Fermi sea term (M a ) this requires the use of Eq. IA3I 



It is useful to define 



t*ni J n ai ni,nj(J > ) = / dr u(r - z x P)$*, „, (r )$ ni ,„ 2 (r) 

(A9) 

« ni ,n. ; ni ,n' 2 (P) = fl(P)*^ in4 (« X P)$ ni ,„ 2 (z X P), 

(A10) 



where u(k) = 2n(3/k is the Fourier transform of the 
Coulomb potential. The quantities —u and v are, in fact, 
the direct and exchange matrix elements for the 2DEG. 
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The graphene matrix elements are then 



For the direct contribution, 



1VI X,2;X' ,2' 



M x 

m X,2;X',2' 
TOl, 2;1 /, 2 '(-P) 



M 



1,2:1' 2' 



<5l,l'#2,2'(£l - £2) 

Sp,P'Ssiy 1 ^h,t' L ^ S 2,s'Jt 2 ,t' 2 rni t2 . 1 , i2 ,(P) 

— l [\/2] + ,° +Sn 2 .° + S n' 2 ,0 

1 

^ ^(ni)Sj/(«2)Sft(ni)Si/(Tl2) 
MA fl (n 1 ),A„(n 2 ):A M (n' 1 ),A„(«^)(-P) 

1 

s / ,(ni)s Jli (n 2 )s I/ (n' 1 )s I/ (n 2 ) 

W A fl (7u),A M (772)^(74), A^ri^-P) 

[^ri 2 ,n! ! , 5ii,t3^si,S3 TO Tll!rl ' 1 ; rl 3 in 3 (0) 
~^ni, n'j <5t 2 , t 3 ^2, S3 TO ri 2 ,n(,, 713,713 (0)], 



(0) 



1,0:1,0 



l 0,-l;0,-l 



" 2 t U 0. 0:0,0 + u l,0;l,o] 

-^J\e- x [{Q + P 2 )I {x)-P 2 h{x)] 



E, 



(i) 



1, 



E. 



(2) 



*2, 1:2,1 - — TP2,1;2,1 
+ M 1,1;1,1 + u 2, 0:2,0 + u l,0;l,o] 

- Ay| e -»[(66 + 13P 2 + 2P 4 + P 6 )/o(^) 
-P 2 (23 + 4P 2 + P 4 )/ 1 (x)] 



1 , 



'3,2:3,2 



[ u 3,2:3,2 



+"2,2;2,2 + «3,1;3,1 + u 2,l;2,l] 

13 ' - e~ x [(1398 + 237P 2 + 10P 4 



3072 V 2 

+43P 6 - 6P 8 + P 10 )I (x) 
-P 2 (497 + 84P 2 + 37P 4 - 4P 6 + P s )h(x)} 



h 



-"0,0;0,0 = -P\ 7T e ^oO) 



1 



'1,0:0,-1 



: M 1,0;1,0 



(x)-(2 + P 2 )I 1 (x)], 



where Xo(n) = 
Si(n) — sgn(n). 



|n|, Ai(n) = |n| - 1, s (n) 



1, and 



The self-energy term has an interesting feature in 
graphene. Note that this term involves the direct ma- 
trix elements evaluated at P = 0. In the 2DEG, the 
conservation of internal angular momentum of the exci- 
ton requires that n i. na „ 3 (0) vanish unless n\ = n\. 

Thus M represents a self-energy correction associated 
with single particle states which depends only on the Lan- 
dau level index. In graphene, the analogous conservation 
law requires only that \n\\ = \n' 1 \. The exchange self- 
energy is thus associated both with a Landau level and its 
electron-hole conjugate, with possible off-diagonal terms. 
Such off-diagonal terms contribute to M only if the chem- 
ical potential is sufficiently low (high) and Landau-level 
mixing is sufficiently strong that both Landau levels —n 
and n > 1 are available for the excited electron (hole). 
This circumstance never arises in our calculations, since 
we restrict our attention to L < 1. 



where x — P 2 /4:, and Iq and I\ are the usual modified 
Bessel functions. 

The exchange matrix elements are evaluated for mo- 
mentum P = —Py: 



E 



(0) 



E. 



(i) 



E. 



(2) 



m l,0;l,0 — m 0, -1:0,-1 — m l,0;0,-l 

= ^ w ,o = fPe- 2 * 

TO 2,1;2,1 = ^[«2,1;2,1 + 2«l, ;2,l + «l,0;l,o] 

^P e - 2a: [24 + 16V2 - 4(2 + V2)P 2 + P 4 ] 

TO 3,2:3,2 = j[«3,2;3,2 + 2«2,1;3,2 + «2,l;2,l] 



1536 



Pe- 2x [192(5 + 2V6) - 96(8 + 3\/6)P 2 



+8(27 + 8\/6)P 4 - 4(6 + \/6)P 6 + P 8 ]. 



The interacting contribution AE"(P) to the exciton 
dispersions computed in this article are given in Table HT1 
in terms of certain matrix elements, which we list below. 



The exchange self-energy constants diverge logarithmi- 
cally in the number of filled levels. We introduce a cutoff 
— rir for the index of the lowest filled Landau level. We 
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define the constants 



,(!) 



,(2) = 



E Ki;«,n(0) - < .„ >n (0)] 

n——n c 
1 

E [ m 2,2;„,„(0) - <1;„,„(0)] 
— — n c 
2 

E K3;n,n(0)-4 



n— — n 
2 



(0)] 



n— — n c 

-1 



E Ko ; «,„(0) ~ rn d _ l! _ 1 . ntn (0)} 

n——n c 

-1 

;n,n (0) 



n— — n c 
-1 



n— — n c 




d = 



E <0;n,JO)- E m -l,-l;n,n(°) 
n——n c n— — n c 

-1 

E TO 0,0;n,n(°) - E m 0,0;«,n(0) 



The summation need only be carried out for a, 
cS- 2 \ and b since 



= -<o,o,o(°)=/?\/£ 



:-o = d- 6= -mf )1;0)0 (0) = 



(All) 



Furthermore, we find that a— b — > in the limit n c — > oo. 
This occurs because a and 6 are exchange constants for 
two many-body states connected by the simultaneous re- 
versal of both charge and energy, and particle-hole sym- 
metry is obtained in the limit n c ~ > oo. A finite cutoff re- 
sults in a difference a — b which is nonvanishing but small 
for n c 3> 1. However, recall that in the tight-binding ap- 
proximation, there are both upper and lower cutoffs for 
the kinetic energy due to the discrete nature of the lat- 
tice. The upper cutoff, and thus particle-hole symmetry, 
is lost in passing to the continuum approximation. In 



this calculation, we assume particle-hole symmetry, and 
consequently a — b. We handle the cutoff by using the 
asymptotic form as n c — > oo: 



= b 



4V2 



-\/7r + 7 + ln(n c ) 



(A12) 



E 



i 1 / 2 T{n- 1/2) 1 



where 7 is the Euler constant. 

The cutoff is determined by the fact that, due to the 
4 spin and pseudospin sublevels, there are approximately 
4n c electrons per quantum of magnetic flux. For a mag- 
netic field of 20 T we find n c m 1870, giving a = b « 2.09, 



d « 2.41, a M 



1(2) 



filling 


branch 


dispersion 


70 = 1,2,3,4; L = 1 


AiJ 1 


a + E^+4E^ 




AE 2 


a + Ei 0) 




AE S 


c + E^+h 




AE 4 


c + E^-h 


71 = 4, l z = 1 


AE 1 


aW+E^+iE^ 




AE 2 


+ EP 


72 = 4, l z = 1 


AE 1 






AE 2 


a^+E (2) 


70 < 4; l z = 


AE° 


e + E 4 



TABLE II: Expressions for various branches of the energy 
shift at different filling factors. 



these values are rather insensitive to value of n c due to 
the logarithmic behavior in Eq. IA12I 

Expressions for the various branches of the energy shift 
are given in Table[n] In the absence of particle-hole sym- 
metry (as can arise from next-nearest neighbor hopping,) 
a — b = c — d =/= 0, and the expressions change somewhat; 
the main effect of this small difference will be a partial 
lifting of the degeneracy of AE 2 in the cases 70 < 4, split- 
ting it into three branches when 70 = 2 and two branches 
when 70 = 1, 3. 
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